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Early brain development is a complex and rapid process, the disturbance of which may cause the onset of brain 
disorders. Based on longitudinal imaging data acquired from 6 to 16 months postnatal, we describe a systematic 
trajectory of monkey brain development during late infancy, and demonstrate the influence of phencyclidine 
(PCP) on this trajectory. Although the general developmental trajectory of the monkey brain was close to that 


ae of the human brain, the development in monkeys was faster and regionally specific. Gray matter volume 
Neurodevelopment began to decrease during late infancy in monkeys, much earlier than in humans in whom it occurs in adolescence. 
Primates Additionally, the decrease of gray matter volume in higher-order association regions (the frontal, parietal and 
Template temporal lobes) occurred later than in regions for primary functions (the occipital lobe and cerebellum). White 


matter volume displayed an increasing trend in most brain regions, but not in the occipital lobe, which had a sta- 
ble volume. In addition, based on diffusion tensor imaging, we found an increase in fractional anisotropy and a 
decrease in diffusivity, which may be associated with myelination and axonal changes in white matter tracts. 
Meanwhile, we tested the influence of 14-day PCP treatment on the developmental trajectories. Such treatment 
tended to accelerated brain maturation during late infancy, although not statistically significant. These findings 
provide comparative information for the understanding of primate brain maturation and neurodevelopmental 
disorders. 


Phencyclidine 


© 2014 Elsevier Inc. All rights reserved. 


Introduction 


Understanding brain development during infancy, a critical stage 
with rapid brain changes and vulnerabilities to exogenous influences, 
is essential for the research of pathological mechanisms of neuro- 
developmental disorders. Recent advances in brain imaging techniques 
enable investigating brain maturation in-vivo. However, partly due to 
the long life span of humans, most human developmental studies 
are based on cross-sectional or mixed-longitudinal data that provide 
less accurate observations than longitudinal studies, especially when 
there is a large individual variability in brain development (Casey 
et al., 2005). Moreover, it is impractical to interfere with normal brain 
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development of humans, which impedes the understanding of neuro- 
developmental disorders. 

The rhesus monkey, a widely-used experimental animal with phylo- 
genetic closeness to humans (Bizon and Woods, 2009), has a shorter life 
span than humans. Thus it is a good candidate for longitudinal studies of 
brain development and neurodevelopmental disorders. However, there 
are few comparative imaging studies on monkey brain development, 
and most of the existing studies are cross-sectional (Chen et al., 2013; 
Shi et al., 2013). A pioneering longitudinal imaging study on monkey 
brain development was performed 8 years ago (Malkova et al., 2006). 
This study investigated the growth of total brain volume and total 
white matter volume from birth to late puberty. Due to the limitation 
of the data (based on 1.5 T MRI scanners), the study did not provide im- 
portant information of gray matter and regional brain growth that are 
essential for the characterization of the trajectories of brain 
development. 


66 C. Liu et al. / Neurolmage 107 (2015) 65-75 


On the other hand, when using monkeys as experimental animals, 
we Can use various ways to disturb brain development. Phencyclidine 
(PCP), a noncompetitive N-methyl-D-aspartate (NMDA) receptor an- 
tagonist, is widely used to build animal models of schizophrenia 
(Jones et al., 2011; Mouri et al., 2007). Recently, schizophrenia has 
been regarded as a neurodevelopmental disorder, whose pathology 
may be contributed by disturbances during prenatal and postnatal 
brain maturation (Jaaro-Peled et al., 2009; Lewis and Levitt, 2002; 
Rapoport et al., 2005). Based on the hypothesis, previous studies have 
examined the effects of PCP on brain development in rodents, and 
they suggest that early postnatal PCP treatment is a novel promising 
animal model for schizophrenia (Broberg et al., 2008, 2013; Kjaerby 
et al., 2014; Rasmussen et al., 2007; Wang et al., 2001). However, spe- 
cies differences between rodents and humans constitute a major gap 
for “translational psychiatry”, which monkeys could help to fill. 

Our study focuses on a comprehensive longitudinal analysis of mon- 
key brain development during late infancy, and on the influence of PCP 
on the development by using non-invasive brain imaging techniques 
and linear mixed-effect models. To our knowledge, this is the first longi- 
tudinal MRI-based study that characterizes regional growth of gray and 
white matter, and the first longitudinal DTI-based study that character- 
izes white matter maturation during early brain development of mon- 
keys, as well as the first work to investigate the influence of PCP on 
these developmental trajectories. 


Materials and methods 
Monkeys 


A total of 14 male rhesus macaques (Macaca mulatta), from the 
breeding colonies at the Primate Center of the Kunming Institute of 
Zoology Chinese Academy of Sciences, were scanned at 5 different 
time-points: a baseline scan (at the age of 196 + 8 days postnatal) 
and about 40, 68, 105, and 303 days later (Fig. 1A, details shown in 
Table S1). These scans covered the period of late infancy, which spans 
from ~0.44 years postnatal (the eruption of the last deciduous tooth) 
to ~1.35 years postnatal (the eruption of the first permanent tooth) 
(Smith et al., 1994). For convenience, these scans were named 6 months 
(S6), 7 months (S7), 8 months (S8), 10 months (S10), and 16 months 
(S16), respectively. In S6, S10 and S16, all 14 monkeys were scanned; 
and in S7 and S8, only 10 monkeys were scanned. After the first scan, 
10 monkeys had a 14-day phencyclidine (PCP, Chemsky (ShangHai) 
International Co., Ltd. #20120516) treatment (twice per day). Of these 
10 monkeys, 5 received PCP with a dose 0.15 mg/kg body weight 
(PCP1), and the other 5 monkeys with a dose of 0.3 mg/kg body weight 
(PCP2). The remaining 4 monkeys were controls treated with saline 
solution (0.1 ml/kg, Control). After birth, the monkeys were reared 
with their mothers in the breeding center. At the age of 6 months, 
they were moved into other cages for experiments (two monkeys per 
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cage), and were housed under standard conditions (a 12-hr light/dark 
cycle with light on from 07:00 to 19:00; humidity at 60%, temperature 
at 21 + 2 °C, 5 times/day deliveries of food and milk). This is the stan- 
dard breeding protocol used in our primate center to reduce the adverse 
influence of maternal separation and the risk of naturally occurring ma- 
ternal neglect and abuse. The experiments were conducted in accor- 
dance with the guidelines for the National Care and Use of Animals 
approved by the National Animal Research Authority of China. Detailed 
scan ages and other information for each monkey are presented in 
Table S1. 


Data acquisition 


MRI and DTI data were acquired with a PHILIPS Achieva 3.0 T MR 
unit of The First Affiliated Hospital of Kunming Medical University, 
Kunming, Yunnan, China. Before the scan, monkeys were intramuscu- 
larly injected with atropine (0.05 mg/kg, ShangHai Harvest Pharmaceu- 
tical Co., LTD.). Then the monkeys were deeply anesthetized with 
ketamine (7.5 mg/kg, im., ShengYang Animal Pharmaceutical factory) 
followed by pentobarbital sodium (5 mg/kg for 1.3-1.9 kg body weight 
and 15 mg/kg for 2.2-3.5 kg body weight, XiTang Biological Company). 
The anesthesia was stable for at least 1 h, and the scan duration was less 
than 1 h. Thus, no complementary anesthesia was needed during the 
scan. A three-dimensional sagittal T1-weighted (T1-w) anatomical 
image was acquired using an inversion-recovery prepared 3-D spoiled 
gradient echo (SPGR) pulse sequence (TR = 14 ms; TE = 7 ms; flip 
angle = 8°; field of view = 160 mm; matrix = 320x320; number of 
averages = 2; final voxel resolution = 0.6 x 0.5 x 0.5 mm’; scan dura- 
tion = 8:30 min). DTI data were acquired using a two-dimensional 
echo-planar diffusion-weighted spin-echo sequence (TR = 8 s; TE = 
78 ms; flip angle = 90°; field of view = 120 mm; matrix = 96 x 96; 
slice thickness = 1.25 mm, no gap; number of averages = 3; final 
voxel resolution = 1.25 x 1.25 x 1.25 mm?; b = 1000 s/mm?; 32 direc- 
tions with one non-diffusion weighted image; scan duration = 15 min). 


Template construction 


MRI (T1-w) template construction 

Brain tissue was manually delineated from non-brain tissue by 
trained researchers for each dataset. Field inhomogeneity was bias- 
corrected using FAST in FSL (Smith et al., 2004; Zhang et al., 2001), 
which also provided segmentation maps including tissue class prob- 
ability images and partial volume images. The MRI template was con- 
structed using the Advanced Normalization Tools (ANTs) software 
(Avants et al., 2008; Klein et al., 2009). First, all images were rotated 
to the same orientation by a rigid alignment, demeaned by dividing 
their mean values, and then averaged to create an initial template 
with 0.5 mm? isotropic resolution. The template was iteratively opti- 
mized by transforming all input images to the current template and 
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Fig. 1. Age and number of monkeys for each scan (A); and procedure of template creation (B). Starting at the second scan, dates are represented as the number of days after the first scan. 
The scans are abbreviated as S6 (6 months), S7 (7 months), S8 (8 months), S10 (10 months), and S16 (16 months). S6, S10 and S16 include all 14 monkeys and are used for linear mixed- 


model analysis. Detailed scan and age information of each monkey are presented in Table S1. 
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averaged to create a refined template. The procedure was performed 
using rigid, affine, and nonlinear alignments, and repeated until the av- 
erage image converged for each type of alignment (Fig. 1B). The trans- 
forms from different alignments were concatenated so that the input 
images were only re-sampled once to minimize interpolation errors. 
After the final template was constructed, probabilistic tissue segment 
maps of each image were warped into the template space, smoothed 
with a 1-mm FWHM Gaussian smoothing kernel, and averaged to create 
tissue probability maps. Furthermore, brain parcellation maps were cre- 
ated using the following steps: 1) The template was transformed to the 
T1-w image of the UNC Primate Brain Atlas (Styner et al., 2007) by an 
affine alignment and then a nonlinear alignment. 2) The template was 
labeled by inverse-transforming UNC parcellation maps to the template 
space based on the transformations estimated from the previous align- 
ments. 3) The template label was carefully checked and manually ad- 
justed. The parcellation maps were used for ROI-based tissue volume 
analysis. 


DTI template construction 

Artifacts from eddy currents and motion in the diffusion data were 
corrected using FDT of FSL (Smith et al., 2004). For each dataset, the 
brain was manually delineated from non-brain tissue on the non- 
diffusion weighted image by trained researchers, and then the diffusion 
tensor was reconstructed using FSL. The DTI template was constructed 
using DTI-TK, a spatial normalization toolkit optimized for examining 
white matter morphometry based on DTI data (Zhang et al., 2007). 
DTI-TK uses higher order information of diffusion tensor images for 
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spatial normalization and has been used to generate high-quality DTI 
templates of adult humans (Zhang et al., 2011), human infants (Wang 
et al., 2011), and adult primates (Adluru et al., 2012). The procedures 
and parameters used to generate a monkey DTI template were similar 
to those employed in a previous study (Adluru et al., 2012). All input 
tensor images were rotated to the same orientation by a rigid alignment, 
and then averaged using a Log-Euclidean mean to create an initial tem- 
plate with a 0.5-mm:? isotropic resolution. Log-Euclidean tensor averag- 
ing preserved white matter orientation with minimal blurring. The 
template was iteratively optimized using rigid and then affine align- 
ments, and repeated until the average image converged for each type 
of alignment. Finally, the template was iteratively optimized using non- 
linear alignments. 


Brain development analysis 


Total brain volume 

Total brain volume was estimated for each subject based on the 
manually de-skulled structural images, and analyzed using linear 
mixed-effect models (see below). 


ROI-based tissue volume analysis 

To avoid bias to particular time points, all structural images were 
transformed to the T1-w template by an affine alignment and then by 
a non-linear alignment. Gray matter and white matter partial volume 
images for each dataset were warped into the template space using 
the transformations from the previous alignments, and modulated by 
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Fig. 2. MRI-based template. The first three columns illustrate the template with corresponding tissue probability maps and parcellation maps. The fourth column shows the representative 


images of one individual monkey. 
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multiplying the Jacobian determinants of the deformation field to com- 
pensate for the expansion/reduction from the nonlinear transformation 
(Good et al., 2001). Finally, tissue volumes for gray and white matter in 
each brain region were estimated based on wrapped partial volume im- 
ages. Besides, the ROI-based analysis was also performed in the native 
space by inverse-transforming the labels of brain regions into the im- 
ages of each monkey to avoid the interpolation and smoothing. The re- 
sults in the native space were similar to the results in the template 
space. 


Voxel-based morphometry (VBM) 

Similar to the ROI-based analysis, gray matter and white matter par- 
tial volume images were warped into the template space and modulat- 
ed by multiplying the Jacobian determinants of the deformation field. 
Wrapped partial volume images were smoothed with a 1-mm FWHM 
Gaussian smoothing kernel and used for the linear mixed-effect analy- 
sis. The analysis was restricted to a gray matter/white matter mask 
from the MRI template with corresponding tissue probabilities greater 
than 20% to exclude voxels with low tissue probabilities (mainly caused 
by spatial smoothing). Finally, the pre-processed images were used for 
linear mixed-effect analysis. Because the ROI-based volume analysis 
provided information similar to that by VBM, the results of VBM were 
complementary and presented in Fig. S1. 


Template 


Deformation-based morphometry (DBM) 

After the affine alignment accounting for the global brain size 
and shape difference, subject brain images of each time point were 
nonlinearly transformed to corresponding images of other time points 
using ANTs (Avants et al., 2008). The Jacobian determinants, which rep- 
resented the local shape difference at each voxel between two time 
points, were then computed for the statistical analysis. 


Tract-based spatial statistics (TBSS) 

TBSS was an observer-independent method for group-wise compar- 
isons of DTI-derived measures (Smith et al., 2006). TBSS comprised 
three steps: 1) normalizing input images to the template space; 2) creat- 
ing a mean fractional anisotropy (FA) skeleton and projecting data to 
the skeleton; 3) statistical analysis. The first step was performed using 
DTI-TK (tensor-based spatial normalization) to register DTI images to 
the template space, because DTI-TK provided better normalization for 
DTI images than the FA-based method implemented in FSL (Bach 
et al., 2014; Zhang et al., 2007). The second step was performed using 
the standard TBSS protocol in FSL (Smith et al., 2004, 2006). FA images 
were calculated from spatial-normalized tensor images and averaged 
to create a mean FA image as well as a mean FA skeleton under an FA 
threshold over 0.2. A lower threshold (FA = 0.15) was also tested due 
to the low FA of infant monkeys, and the results were similar to those 


Individual 


Fig. 3. DTI-based template. The first three columns show the properties of the final template, and the fourth column shows the properties of an individual monkey image. Fractional an- 
isotropy (FA), trace, axial diffusivity (AD), radial diffusivity (RD) and principal diffusion direction images (PD) are presented in each row. The units of trace, AD, and RD are mm?/s. The 
inferior cingulum mask required for tract-based spatial statistics (TBSS) is also shown on the sagittal slice of the FA image (in red color). 
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under the threshold of 0.2. An inferior cingulum mask was manually 
drawn on the template space to facilitate FA projection because the in- 
ferior cingulum tract had no well-defined search direction for the pro- 
jection due to its tubular skeleton topology (Smith et al., 2006). After 
the projection, linear mixed-effect models were used to analyze the 
skeletonized FA images. Concurrently, TBSS of other diffusion-derived 
measures, including mean diffusivity (MD), axial diffusivity (AD) and 
radial diffusivity (RD) were also performed using the same approach. 


Statistical analysis 

For the total brain volume and ROI-based tissue volume, time- 
related change and PCP-related change were analyzed using linear 
mixed-effect models by the NLME package of the R language (Pinheiro 
et al., 2007). Main (fixed) effects of the time, the group factor (control, 
PCP1 and PCP2) and the interaction between the time and the group 
factors (time-group interaction) were included in the analysis; and 
the random effect (subject-specific intercept) of each monkey was 
also modeled to account for repeated measurements. Time could be 
modeled as a categorical factor to test the difference between different 
time points (M1), or it could be modeled as a continuous variable to 
test the linear/quadratic trend (M2), resulting in two similar but differ- 
ent models. For the M2, the quadratic term of time was not significant in 
any of the measures; thus, a linear fit was calculated. Besides, more 
complex structures of variance, correlation and random effect were 
also modeled, including heterogeneous variance, correlation structure 
with an autoregressive process of order 1, and random effect with 
both subject-specific intercept and slope. These models were compared 
by likelihood-ratio tests and information criteria (AIC and BIC), but none 
of these complex models outperformed the mixed-effect models M1 
and M2. Thus, M1 and M2 were chosen for the final mixed-effect 
model analysis, which were followed by Tukey post-hoc tests for any 
significant factors. 

Similar models (M1 and M2) were fitted for VBM and TBSS by 
SurfStat (Worsley, 2009), which was capable of fitting mass-univariate/ 
multivariate linear mixed-effect models for the volumetric data. As the 
F tests for mixed effect models were not featured in SurfStat, different 
contrasts were designed to test different factors. 

DBM was calculated from a nonlinear alignment between the 
images of two scans. For each pair of time points, a voxelwise one- 
sample permutation test was performed on the maps of Jacobian deter- 
minants to determine if there was a significant expansion or reduction 
of the brain using the randomise of FSL (Smith et al., 2004). 


Results 
Templates of infant monkeys 


To avoid any bias to a particular time point and to improve the qual- 
ity of templates, we created high-quality MRI and DTI templates using 
the data of all scans. The templates were a pre-requisite for the subse- 
quent brain imaging analysis, and served as a target space to which in- 
dividual brains were spatially normalized. Both the MRI-based template 
and the DTI-based template are provided in the Supplementary Data for 
downloading and usage. 


MRI-based template 

As illustrated in Fig. 2, the template had a superior contrast between 
different tissues and a better signal-to-noise ratio when compared to 
individual brains. The greater contrast and improved delineation of 
anatomical structures in the template demonstrated a good alignment 
of individual monkey brains. Moreover, the template contained brain 
parcellation maps derived from the UNC Primate Brain Atlas. Thus, we 
can use the template to label each brain region and to conduct ROI- 
based analysis. 


DTI-based template 

The properties of the DTI-based template and a representative indi- 
vidual example, including fractional anisotropy (FA), trace (trace = 3x 
mean diffusivity), axial diffusivity (AD), radial diffusivity (RD), and prin- 
cipal diffusion direction images (PD) are shown in Fig. 3. We also drew 
manually an inferior cingulum mask on the template space to facilitate 
its skeleton-projection for the following TBSS analysis. Like the MRI- 
based template, the DTI-based template had a better contrast than indi- 
vidual brains, which demonstrated a good alignment. 


Brain development and the effect of PCP 


We characterized the rapid development of the monkey brain and 
the influence of PCP during late infancy based on longitudinal brain im- 
aging data. The results across S6, S10 and S16 are described in this sec- 
tion, as 1) the three time points showed evident changes in brain 
maturation, 2) S10 and S16 were far apart from the baseline scan (S6) 
and suffered less from the variability of the ages of different monkeys, 
3) the three time points had all 14 monkeys scanned, which allowed 
more solid statistical inference. 


Total brain volume 

We calculated the total brain volume based on individual de- 
skulled structural images. Mixed-effect model analysis showed a sig- 
nificant growth of brain volume during late infancy (F(2,22) = 
154.63, p < 0.001): on the average, the 10-month-old brain and 
the 16-month-old brain were 1.062 and 1.092 x larger than the 6- 
month-old brain, respectively. However, the rate of brain growth grad- 
ually slowed down: an increase of 6.21% occurred from 6 to 10 months 
(adjusted p < 0.001), and was then followed by another increase of 
2.81% from 10 to 16 months (adjusted p < 0.001) (Fig. 4). 

Although the PCP-treated groups had higher growth rates (11.15% in 
PCP1, 8.86% in PCP2) than the control group (7.13%) from 6 to 16 
months (Fig. 4), neither the group factor (F(2,11) = 0.38, p = 0.69) 
nor the time-group interaction (F(4,22) = 2.58, p = 0.07) had signifi- 
cant effects on the changes of total brain volume. 
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Fig. 4. Total brain volume changes. The percent changes of total brain volume of each scan 
against the baseline scan (S6) are presented in mean + SEM for each group. Mixed-effect 
model analysis (based on absolute volume) shows significant time-related changes 
(F(2,22) = 154.63, p < 0.001), but non-significant effects of the group factor (F(2,11) = 
0.38, p = 0.69) and non-significant time-group interaction (F(4,22) = 2.58, p = 0.07). 
Asterisks indicate the difference of volume between the scans based on post-hoc tests: 
*p < 0.05, **p < 0.01, and ***p < 0,001. Detailed statistics of mixed-model analyses are pre- 
sented in Table S2. 
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ROI-based tissue volume analysis 

Development of the monkey brain during late infancy was charac- 
terized by reduced gray matter volume and increased white matter vol- 
ume (Figs. 5 and 6). However, the developmental trend of tissue volume 
was regionally specific. Gray matter volume in the occipital lobe, cere- 
bellum, and insular cortex decreased significantly after S6. Gray matter 
in other cortices, including the frontal, parietal, and temporal lobes, and 
cingulate cortex, decreased significantly only after S10. White matter 
volume displayed an overall increasing trend in most brain areas, but 
not in the occipital lobe (F(2,22) = 1, p = 0.38). A non-significant in- 
crease of white matter volume was also found in the parietal lobe 
(F(2,22) = 2.52, p = 0.10). The increase of white matter volume from 
S6 to S10 was significant in the cerebellum and corpus callosum, but 
not significant in the frontal lobe, temporal lobe, and cingulate cortex. 
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Although neither the group factor nor the time-group interaction 
had significant effects on any of the above changes, the PCP groups 
displayed a non-significant acceleration of the decrease in gray matter 
volume and of the increase in white matter volume after S10 in most 
brain regions, including the frontal, parietal, and temporal lobes, cere- 
bellum, and cingulate cortex (Figs. 5 and 6). All detailed statistics of 
mixed-model analyses and the absolute volume of the above measures 
are presented in Table S2. 


Voxel-based morphometry (VBM) and Deformation-based morphometry 
(DBM) 

ROI-based analysis provided measures of large areas, whereas VBM 
enabled us to compare changes of tissue volume at the voxel level. Con- 
sistent with ROI-based analysis, VBM revealed reduced gray matter 
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Fig. 5. ROI-based gray matter volume analysis. Percent changes of gray matter volume of each scan against the baseline scan (S6) are presented in mean + SEM for each group. Mixed- 
effect model analysis (based on absolute volume) shows significant time-related changes in all brain regions. Neither the group factor nor the time-group interaction is significant in any 
brain regions. Asterisks indicate differences of volume between the scans based on post-hoc tests: *p < 0.05, **p < 0.01, and ***p < 0.001. Detailed statistics of mixed-model analyses are 


presented in Table S2. 
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Fig. 6. ROI-based white matter volume analysis. Percent changes of white matter volume of each scan against the baseline scan (S6) are presented in mean + SEM for each group. Mixed- 
effect model analysis (based on absolute volume) shows significant time-related changes in most brain regions, but not in the occipital and parietal lobes. Neither the group factor nor the 
time-group interaction is significant in any brain region. Asterisks indicate the difference of volume between the scans based on post-hoc tests: *p < 0.05, **p < 0.01, and ***p < 0.001. De- 


tailed statistics of mixed-model analyses are presented in Table S2. 


volume and increased white matter volume across the brain during late 
infancy (Fig. S1). VBM also showed regional variations: white matter 
changes in the occipital and parietal lobes were less obvious than the 
changes in other brain regions. Neither the group factor nor the time- 
group interaction had significant effect in the changes of VBM. 

DBM was an effective technique to detect morphological differences 
based on the deformation field derived from nonlinear alignments. Con- 
sistent with the above findings, more expansion was found in the white 
matter than the gray matter, whereas more reduction was found in the 
gray matter than the white matter (Fig. 7). 


Tract-Based Spatial Statistics (TBSS) 
Besides volume and shape changes, we characterized changes 
in white matter microstructure based on DTI. TBSS showed increased 


fractional anisotropy (FA) and decreased diffusivity (MD, RD and AD) 
across all major tracts of the brain during late infancy (Fig. 8). The 
decrease in RD (diffusivity perpendicular to the axon) was more evident 
than that of AD (diffusivity along the axon). Thus, changes in FA and MD 
were driven mainly by the reduction in RD. All of these changes were 
more obvious at S16 than S10, suggesting continuous development of 
the white matter microstructure during late infancy. None of the mea- 
sures showed any significant effect of PCP on brain development. 


Discussion 
This study demonstrates that 1) the volume of gray matter decreases 


in monkeys during late infancy, much earlier than in humans; 2) the 
volume of white matter increases and is associated with an increase in 
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Deformation-based morphometry 
A: S10 - S6 B: S16 - S10 C: S16 - S6 


P < 0.05, FWE and TFCE corrected 


Fig. 7. Deformation-based morphometry. The red color represents significant expansion and the blue color represents significant reduction based on the permutation test implemented in 
the randomise of FSL (Smith et al., 2004). 
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P < 0.05, FDR corrected 


Fig. 8. Tract-Based Spatial Statistics (TBSS). Fractional anisotropy, mean diffusivity, axial diffusivity, and radial diffusivity are analyzed by mass-univariate linear mixed-effect models. The 
first three figures of each property represent the time-related contrasts of M1 where the time is modeled as categorical factors; the fourth figure represents the time-related contrasts of 
M2 where the time is modeled as a continuous variable; and the fifth figure represents all contrasts related to the group factor and the time-group interaction. All skeletonized results 
(yellow and azure) are filled into local tracts (red and blue) by the tbss_fill of FSL for a better visualization (Smith et al., 2006). The red and yellow colors represent significant increases, 
and the blue and azure colors represent significant decreases in these properties. 
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Fig. 9. Summary of human and monkey brain development from early infancy to puberty (adolescence). Adopted from a recent study (Sakai et al., 2011), these developmental stages are 
defined based on combined dental eruption and sex maturation for inter-specific comparisons (Plant, 1994, 2001; Plant and Barker Gibb, 2004; Smith et al., 1994; Terasawa and Fernandez, 
2001). The time scale of monkeys is re-scaled to match that of humans for a better comparison. These developmental trajectories are depicted mainly based on longitudinal and mixed- 
longitudinal data (Brouwer et al., 2012; Giedd, 2004; Giedd et al., 1999; Lebel and Beaulieu, 2011; Malkova et al., 2006) with the cross-sectional data as supplements (Geng et al., 2012; 
Knickmeyer et al., 2008, 2010; Lenroot et al., 2007; Shi et al., 2013). The phases based on (mixed-)longitudinal data are marked with solid outlines, the cross-sectional data with dashed 
outlines, and our longitudinal data with red solid outlines. The color bar from white to black represents an increase in measures, and vice versa. 
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fractional anisotropy and a decrease in diffusivity; and 3) the matura- 
tion of gray and white matter are regionally specific; and 4) PCP tends 
to accelerate maturation after 10 months of age, although not statistical- 
ly significant. 

We observed a rapid growth of monkey brain volume during early 
development, consistent with a pervious study (Malkova et al., 2006). 
However, the growth rate began to slow down during late infancy, 
as 6.21% of the expansion occurred between 6 and 10 months and 
only 2.81% between 10 and 16 months (Fig. 4). The fast expansion is 
mainly caused by the significant increase of white matter volume 
(Fig. 6), while the later decline in the expansion rate can be partly 
explained by a significant decrease of gray matter volume after 10 
months (Fig. 5). Compared with monkey brain development, human 
brain volume increases 100% during the first year of life, and then 
keeps cubic growth into adolescence (Fig. 9) (Giedd, 2004; Giedd 
et al., 1999; Gilmore et al., 2007), which indicates that humans have a 
less mature brain at birth and a longer period of early rapid brain 
growth than monkeys (Gibson and Petersen, 1991). 

The changes in the volume of gray and white matter volume of 
monkey brains were regionally specific (Figs. 5 and 6). The occipital 
lobe and cerebellum (presumably important for primary functions) 
had a significant decrease in gray matter volume before 10 months, 
whereas the frontal, temporal, and parietal lobes (associated with func- 
tional integrations) had significant decreases only after 10 months. As to 
white matter, its volume showed an increasing trend in most brain re- 
gions during late infancy, but not in the occipital lobe. The developmen- 
tal pattern demonstrates that brain regions for primary functions 
mature earlier than brain regions associated with functional integration. 
In humans, brain maturation shows more complex regional variations 
(Fig. 9). Gray matter volume in brain regions for functional integration, 
including the frontal, temporal, and parietal lobes, decreases only after 
adolescence and the change tends to follow a trajectory of an “inverted 
U curve” with peaks at different times, whereas the gray matter of the 
occipital lobe keeps increasing without significant decline or leveling 
(Giedd, 2004; Giedd et al., 1999). Compared with monkeys, the more 
complex and delayed brain development of humans allows a stronger 
influence of postnatal experiences on brain structures, which may con- 
tribute to the emergence of human-specific cognitive abilities and social 
behaviors. 

In addition to volume change, we also characterized the develop- 
ment of white matter microstructure using DTI. We observed an overall 
increase in FA and a decrease in diffusivity (MD, RD, and AD). Previous 
studies revealed increased RD and decreased FA in mice that were ge- 
netically modified or treated with cuprizone to be myelin-deficient 
(Nair et al., 2005; Song et al., 2002, 2005). As demyelination contributes 
to a decrease of FA and an increase of RD, we infer that myelination may 
be the main reason for the change of FA and RD during late infancy. We 
also observed a decrease in AD, which is consistent with previous 
findings in humans (Fig. 9) (Forbes et al., 2002; Mukherjee et al., 
2001; Sadeghi et al., 2013; Schneider et al., 2004; Zhang et al., 2005). 
Contrary to human development, a previous cross-sectional DTI study 
in monkeys showed an increase of AD from late infancy to young adult- 
hood (Shi et al., 2013). The above evidence indicates that the change of 
AD in monkeys may be nonlinear after late infancy, which is worth fur- 
ther longitudinal validation (Fig. 9). 

We also investigated the effect of PCP on these trajectories with 
the hope of using monkeys as an animal model for schizophrenia 
(Jones et al., 2011; Mouri et al., 2007). Despite their trajectories being 
similar to those of the control group, the PCP-treated groups had an ac- 
celerated rate of development after 10 months (Figs. 4-6). However, 
none of the PCP-related changes was significant. Two possible reasons 
may account for the non-significance. The first may be the small sample 
size (14 monkeys). However, a longitudinal study like ours has a higher 
power and requires fewer samples than cross-sectional and mixed- 
longitudinal studies. In addition, we are mainly interested in major 
PCP-treatment effects that are detectable even with a small sample 


size. The second reason may be that the time of PCP treatment was 
beyond the critical phases of the development of the glutamatergic 
circuitry in the prefrontal cortex that PCP primarily affects: one is 
from late gestation to the first 3 months postnatal, when the number 
of excitatory synapses (glutamatergic) substantially increases; the 
other from around 15-18 months to young adulthood, when the num- 
ber of excitatory synapses begins to decline until the stable adult level is 
achieved (Hoftman and Lewis, 2011). In our study, the time points of 
PCP treatment and scans were within the stable period from 3 to 15- 
18 months, which partly explains the non-significant effect of PCP. 
Therefore, it is worthwhile to study if the effect of PCP becomes signifi- 
cant or disappears when these monkeys reach adulthood, and whether 
these changes can influence cognitive performance that are demanding 
for infant monkeys. 

Besides the issue of PCP, we cannot investigate sex differences in 
brain development because the monkeys involved in our study were 
all males. Although the general developmental trajectory of males 
may be applicable to females, gender differences in early brain develop- 
ment may exist and are worth future investigation. Furthermore, our 
longitudinal data only cover development during late infancy, a time 
scale that needs to be expanded in the future. 

Supplementary data to this article can be found online at http://dx. 
doi.org/10.1016/j.neuroimage.2014.11.056. 
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